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1 , Abstract 

O ' 

^ ' We investigate theoretically whether it is feasible to detect 77- and w-nucleus bound 

a : 

states. As well as the closed shell nuclei, 16 0, 40 Ca, 90 Zr and 208 Pb, we also investigate 
^ . 6 He, n B and 26 Mg, which are the final nuclei in the proposed experiment involving the 

(d, 3 He) reaction at GSI. Potentials for the 77 and to mesons in these nuclei are calculated 
in local density approximation, embedding the mesons in the nucleus described by solv- 
ing the mean-field equations of motion in the QMC model. Our results suggest that one 
should expect to find 77- and w-nucleus bound states in all these nuclei. 
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The study of the properties of hadrons in a hot and/or dense nuclear medium is one of the 
most exciting new directions in nuclear physics. In particular, the medium modification of the 
light vector (p, uo and 0) meson masses has been investigated extensively by many authors 
It has been suggested that dilepton production in the nuclear medium formed in relativistic 
heavy ion collisions, can provide a unique tool to measure such modifications as meson mass 
shifts. For example, the experimental data obtained at the CERN/SPS by the CERES [0 
and HELIOS || collaborations has been interpreted as evidence for a downward shift of the p 
meson mass in dense nuclear matter 0]. To draw a more definite conclusion, measurements of 
the dilepton spectrum from vector mesons produced in nuclei are planned at TJNAF || and 
GSI[§ 

Recently, a new, alternative approach to study meson mass shifts in nuclei was suggested 
by Hayano et al. 0. Their suggestion is to use the (d, 3 He) reaction to produce rj and u 
mesons with nearly zero recoil. If the meson feels a large enough, attractive (scalar) force 
inside a nucleus, the meson is expected to form meson-nucleus bound states. Hayano et al. [Q 
have estimated the binding energies for various ?7-mesic nuclei. They have also calculated some 
quantities for the u meson case. However, they used an ^-nucleus optical potential calculated 
to first-order in density, taking as input the 77-nucleon scattering length. In this article, we use 
an alternative, self-consistent method to study whether it is possible to form 77- and cu-nucleus 
bound states in 16 0, 40 Ca, 90 Zr and 208 Pb, as well as 6 He, n B and 26 Mg. The latter three 
nuclei correspond to the proposed experiments at GSI using the (d, 3 He) reaction - i.e., the 
reactions, 7 Li (d, 3 He) ^He, 12 C (d, 3 He) J^B and 27 A1 (d, 3 He) 2 ^Mg. 

In earlier work we addressed the question of whether quarks play an important role in finite 
nuclei |9|, [10], [Tl[] . The quark- meson coupling (QMC) model |12|, §, which is based explicitly 
on quark degrees of freedom, is probably one of the most appropriate models to study whether 
meson-nucleus bound states are possible. The model has been able to describe successfully the 



static properties of both nuclear matter and finite nuclei |10|, [13[], as well as meson properties 



in the nuclear medium [fLlj| . Thus, the model is ideally suited to treat a bound meson and the 
nucleons in a nucleus on the same footing. In this study, we investigate the possible formation 
of the 77- and cj-nucleus bound states due to downward shifts of the masses. We will use QMC- 



I 111 111, where the effective, isoscalar- vector to field, which is off mass-shell and mediates the 



1 



interactions among nucleons, is distinguished from the physical (on mass-shell) u meson which 
is produced inside the nuclei by the above mentioned experiments. 

One of the most attractive features of QMC is that, in practice, it is not significantly 
more complicated than Quantum Hadrodynamics (QHD) |14|], although the quark substructure 
of hadrons is explicitly implemented. Furthermore, it produces a reasonable value for the 
nuclear incompressibility. A detailed description of the Lagrangian density, and the mean-field 
equations of motion needed to describe a finite nucleus, is given in Refs. [0, (TT 
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At position r in a nucleus (the coordinate origin is taken at the center of the nucleus), 
the Dirac equations for the quarks and antiquarks in the rj and to meson bags are given by [[□] : 



i 1 ■ d x - (m q - V a (r)) =F 7° (v u (f) + l -V p {r) 



il-d x - (m q - V a (r}) =F 7° \V u (f) - -V t 




0. 



o. 



0. 



(2) 
(3) 



[ry ■ d x - m 8 ] ip s {x) (or ip s (x)) 

(Note that we have neglected a possible, very slight variation of the scalar and vector mean- 
fields inside the meson bag due to its finite size ||.) The mean-field potentials for a bag 
centered at position r in the nucleus, which will be calculated self-consistently, are defined by, 
v <r{r) = #Mr)>K>(r) = gfa(r) and V p {r) = g q p b(f), with g*,g* and g g p being, respectively, 
the corresponding quark and meson-field coupling constants. Here we assume that the current 
masses are given as m q = m u = = rriu = vng. Furthermore, we have assumed that the a, uj 
and p fields only interact directly with the nonstrange quarks and antiquarks [|11|]. The mean 
meson fields at position r in the nucleus are calculated self-consistently by solving Eqs. (23) - 
(30) of Ref. 0. 

Hereafter we use the notation, ujb, to specify the physical, bound uo meson, in order to 
avoid confusion with the isoscalar-vector uj field appearing in QMC. The static solution for the 
ground state quarks or antiquarks in the rj and uj meson bags may be written as: 



ipf(x) — Nfe %eft l R iipf{x), (for j = r],UB and f — u,u,d,d,s,s), 



(4) 



where Nf and ipf(x) are respectively the normalization factor and corresponding spin and 
spatial part of the wave function [|ll|]. The bag radius in medium, R* (j = i],ujb), which 



depends on the hadron species in which the quarks and antiquarks belong, will be determined 
self-consistently through the stability condition for the (in-medium) mass of the meson against 
the variation of the bag radius. (See Eq. fliCf ) below.) The eigenenergies for the quarks, in units 
of l/Rj, are given by 




tr q (f) ± r* (v u (f) + \v p {f^ 



r s[ r) ' v - 

e s (r) = e s (r) = n s (r), (5) 



where £l*(r) = Jx 2 , + (R*m*) 2 and Q s (r) = y x\ + (R*m s ) 2 with m* = m q — g^a(r)(q = 
u, u, d, d). The bag eigenfrequencies, x q and x s , are determined by the usual, linear boundary 
condition ||. 

Next, we consider the rj and uj p meson masses in the nucleus. The physical states of the 
rj and ujb mesons are the superpositions of the octet and singlet states: 

i = £ 8 cos P y - 6 sin 9 P y, £' = £ 8 sin 9 P y + 6 cos 9 P y, (6) 

with 

£i = —/= (uu + dd + ss), £s = — p {uu+ dd — 2ss), (7) 

where (£, £') denotes (rj, rj') or (</>, lub), with the mixing angles 9 P or 9y, respectively [ [L5[] . 
Then, the masses for the rj and cus mesons in the nucleus at the position r, are self-consistently 
calculated by: 



2\apn* (r) + bpttJr)] - z„ 4 
m;(r) = 1 P 91 J ^ " K )l V - + -7rR* v *B, (Sj 



*3 ; 



V 



2[a^r q (r) + b^U?)} - z«„ , 4 . ; 



dm*(r) 



dRj 
with 



< B <f) = q R : - +i«KLB, (9) 

= 0, U = V,Ub), (10) 

-Rt =R* 



1/2 /2 1 

ap,y = cos 6» P>y - y - sin fl^y, 6 Py = y - cos 9 P y + -j= sin 9 P y. (11) 

In practice, we use 9 P = —10° and 9y = 39° [[HJ, neglecting any possible mass dependence 
and imaginary parts. We also assume that the values of the mixing angles do not change in 



medium, although this is possible and merits further investigation. In Eqs. (H) and @, z v and 
z^ B parameterize the sum of the center-of-mass and gluon fluctuation effects, and are assumed 
to be independent of density 0. 

In this study, we chose m q = 5 MeV and m s = 250 MeV, for the current quark masses, and 
Rn = 0.8 fm for the bag radius of the nucleon in free space. Other inputs, parameters, and some 
of the quantities calculated in the present study, are listed in Table [3]. The coupling constants, 
g%, g% and g q , are adjusted to fit the saturation energy and density of symmetric nuclear matter 
and the bulk symmetry energy. Note that none of the results for nuclear properties depend 
strongly on the choice of the other parameters - for example, the relatively weak dependence 
of the final results on the values of the current quark mass and bag radius is shown explicitly 
in Refs. || p70f| . The parameters at the hadronic level associated with the core nucleus are 
summarized in Table 0. The value of the a mass for finite nuclei is obtained by fitting the 
r.m.s. charge radius of 40 Ca to the experimental value, r c h( 40 Ca) = 3.48 fm |ID[ . For more 
details and explanations of the model parameters, see Refs. || |IOfl . 



Table 1: Physical masses fitted in free space, free space full widths, T, the bag parameters, 
z, and the bag radii in free space, R. The quantities with an asterisk, are those quantities 
calculated at normal nuclear matter density, po — 0.15 fm" 3 . They are obtained with the bag 
constant, B = (170 MeV) 4 , current quark masses, m u = rrid = 5 MeV and m s = 250 MeV. 
Note that the free space width of the rj meson is 1.18 keV fl~5H . 





mass (MeV) 


T (MeV) 


z 


R (fm) 


m* (MeV) 


R* (fm) 


N 


939.0 (input) 




3.295 


0.800 (input) 


754.5 


0.786 


V 


547.5 (input) 


(input) 


3.131 


0.603 


483.9 


0.600 


UJ B 


781.9 (input) 


8.43 (input) 


1.866 


0.753 


658.7 


0.749 



Through Eqs. (P - (|TT|) we self-consistently calculate effective masses, m*(r) and rrf^^f) 
at the position r in the nucleus. Because the vector potentials for the same flavor of quark and 
antiquark cancel each other, the potentials for the i] and ujb mesons are given respectively by 
m*(r) — and m* (r) — m u)B) where they will depend only on the distance from the center 
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Table 2: Parameters at the hadronic level (masses and coupling constants of mesons and photon 
for finite nuclei) |T0| . 



field 


mass (MeV) 


g 2 /4n (e 2 /47r) 


o 


418 


3.12 


UJ 


783 


5.31 


p 


770 


6.93 


A 





1/137.036 



of the nucleus, r = \r\. Before showing the calculated potentials for the i] and ujb, we first 
show in Fig. [I] their effective masses and those calculated within an SU(3) quark model basis, 
uj = -^{uu + dd) (ideal mixing) and t] 8 = in Eq. (|7j), in symmetric nuclear matter. One can 
easily see that the effect of the singlet-octet mixing is negligible for the ujb mass in matter, 
whereas it is important for the rj mass. 

As an example, we show the potentials for the mesons in 26 Mg and 208 Pb in Fig. |2|. Note 
that the actual calculations for 6 He, n B and 26 Mg are performed in the same way as for the 
closed shell nuclei, 16 0, 40 Ca, 90 Zr and 208 Pb. Although 6 He, n B and 26 Mg are not spherical, 
we have neglected the effect of deformation, which is expected to be small and irrelevant for the 
present discussion. (We do not expect that deformation should alter the calculated potentials 
by more than a few MeV near the center of the deformed nucleus, because the baryon (scalar) 
density there is also expected to be more or less the same as that for a spherical nucleus - 
close to normal nuclear matter density.) The depth of the potentials are typically 60 and 130 
MeV for the i] and ujb mesons, respectively, around the center of each nucleus. In addition, we 
show the calculated potentials using QMC-II |Tl| in Fig. ^|, for 208 Pb, in order to estimate the 
ambiguities due to different versions of the QMC model. At the center of 208 Pb, the potential 
calculated using QMC-II is about 20 MeV shallower than that for QMC-I. 

Now we are in a position to calculate single-particle energies for the mesons using the 
potentials calculated in QMC. Because the typical momentum of the bound uj is low, it should 
be a very good approximation to neglect the possible energy difference between the longitudinal 
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and transverse components of the uo |T^|. Then, after imposing the Lorentz condition, <9 M M = 0, 
solving the Proca equation becomes equivalent to solving the Klein-Gordon equation, 



V 2 +K 2 -mf(r) &(f) = 0, ti = V,UB), 



:i2) 



where Ej is the total energy of the meson. An additional complication, which has so far been 
ignored, is the meson absorption in the nucleus, which requires a complex potential. At the 
moment, we have not been able to calculate the imaginary part of the potential (equivalently, 
the in-medium widths of the mesons) self-consistently within the model. In order to make a 
more realistic estimate for the meson-nucleus bound states, we include the widths of the rj and 
ojb mesons in the nucleus by assuming a specific form: 



i r 



m*(r) = m*{r) — - (rrij — m*(r))jj + T 



rrij [r ) 



5 rj(r) 



(13) 
(14) 



where, rrij and Tj are the corresponding masses and widths in free space listed in Table [3], and 
7j are treated as phenomenological parameters to describe the in-medium meson widths, T*(r). 
According to the estimates in Refs. J7L 17], the widths of the mesons in nuclei and at normal 
nuclear matter density are r* ~ 30 — 70 MeV and T* s ~ 30 - 40 MeV 0, respectively. 
Thus, we calculate the single-particle energies for several values of the parameter, 7^, which 
cover the estimated ranges. 

From Table [I] and the calculated density distributions one can obtain the corresponding 
widths at normal nuclear matter density, as well as in the finite nuclei. Because of the recoilless 
condition for meson production in the GSI experiment 0, f§, we may expect that the energy 
dependence of the potentials would not be strong Thus we actually solve the following, 
modified Klein-Gordon equations: 



V 2 + £j-mf(r) 0,(r-) = O, U = V,Ub). 



(15) 



This is carried out in momentum space by the method developed in Ref. [IS]. To confirm the 
calculated results, we also calculated the single-particle energies by solving the Schrodinger 
equation. Calculated single-particle energies for the r\ and ujb mesons, obtained solving the 
Klein-Gordon equation are respectively listed in Tables |3| and [|. 
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Table 3: Calculated rj meson single-particle energies, E = Re{E v — m v ), and full widths, T, 
(both in MeV), in various nuclei, where the complex eigenenergies are, E v = E + m v — iT/2. 
See Eq. ( |l~3"D for the definition of j v . Note that the free space width of the rj is 1.18 keV, which 
corresponds to 7^ = 0. 





7,=0 


7,=0.5 


7,=l-0 




E Y 


E Y 


E Y 

1 j 1 


* 6 Is 
Id 


-33.1 
-8.69 


-32 6 26 7 
-7.72 18.3 


-31 2 53 9 
-5.25 38.2 


^°Ca Is 
lp 
2s 


-46.5 
-27.4 
-6.09 


-46.0 31.7 
-26.8 26.8 
-4.61 17.7 


-44.8 63.6 
-25.2 54.2 
-1.24 38.5 


9°Zr Is 
lp 
2s 


-53.3 
-40.5 
-22.3 


-52.9 33.2 
-40.0 30.5 
-21.7 26.1 


-51.8 66.4 
-38.8 61.2 
-19.9 53.1 


f*Pb Is 
lp 
2s 


-56.6 
-48.7 
-36.3 


-56.3 33.2 
-48.3 31.8 
-35.9 29.6 


-55.3 66.2 
-47.3 63.5 
-34.7 59.5 


°He Is 


-11.4 


-10.7 14.5 


-8.75 29.9 


?B Is 


-25.0 


-24.5 22.8 


-22.9 46.1 


fMg Is 
lp 


-39.2 
-18.5 


-38.8 28.5 
-17.8 23.1 


-37.6 57.3 
-15.9 47.1 
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Table 4: As in Tables |3], but for u meson single-particle energies. In the light of V in Refs. [17], 
the results with 7^ = 0.2 are expected to correspond best with the experiment. 





7^=0 


7^=0-2 


7^=0.4 




E 


r 


E 


r 


E 


r 




Is 


-93.5 


8.14 


-93.4 


30.6 


-93.4 


53.1 




lp 


-64.8 


7.94 


-64.7 


27.8 


-64.6 


47.7 




Is 


-111 


8.22 


-111 


33.1 


-111 


58.1 




lp 


-90.8 


8.07 


-90.8 


31.0 


-90.7 


54.0 




2s 


-65.6 


7.86 


-65.5 


28.9 


-65.4 


49.9 


90 Zr 


Is 


-117 


8.30 


-117 


33.4 


-117 


58.6 




lp 


-105 


8.19 


-105 


32.3 


-105 


56.5 




2s 


-86.4 


8.03 


-86.4 


30.7 


-86.4 


53.4 


208 pu 


Is 


-118 


8.35 


-118 


33.1 


-118 


57.8 




lp 


-111 


8.28 


-111 


32.5 


-111 


56.8 




2s 


-100 


8.17 


-100 


31.7 


-100 


55.3 


£He 


Is 


-55.7 


8.05 


-55.6 


24.7 


-55.4 


41.3 




Is 


-80.8 


8.10 


-80.8 


28.8 


-80.6 


49.5 




Is 


-99.7 


8.21 


-99.7 


31.1 


-99.7 


54.0 




lp 


-78.5 


8.02 


-78.5 


29.4 


-78.4 


50.8 




2s 


-42.9 


7.87 


-42.8 


24.8 


-42.5 


41.9 
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Our results suggest one should expect to find bound rj- and cu-nuclear states as has been 
suggested by Hayano et al. j7|, §]. For the r] single-particle energies, our estimated values lie 
between the results obtained using two different parameter sets in Ref M. From the point of 
view of uncertainties arising from differences between QMC-I and QMC-II, the present results 
for both the single-particle energies and calculated full widths should be no more than 20 % 
smaller in absolute values according to the estimate from the potential for the us in 208 Pb in 
Fig. 0. Nevertheless, for a heavy nucleus and relatively wide range of the in-medium meson 
widths, it seems inevitable that one should find such rj- and u;-nucleus bound states. Note 
that the correction to the real part of the single-particle energies from the width, T, can be 
estimated nonrelativistically, to be of order of ~ T 2 /8m (repulsive), which is a few MeV if we 
use T ~ 100 MeV. 

In future work we would like to include the effect of g-uj mixing, which (within QHD, 



at least) becomes especially important at higher densities JT6|. It will also be important for 
consistency to calculate the in-medium width of the meson within the QMC model and to 
study the energy dependence of the meson-nucleus potential. While the energy dependence of 
the potential felt (for example) by the uo may be quite significant as we move from a virtual u> 
(q 2 ~ 0) to an almost real u> (q 2 ~ m 2 ) [17]], QHD studies in nuclear matter did not reveal a 
strong energy dependence for q 2 near m 2 ][]1| - the region of interest here. Nevertheless, this 
point merits further study in finite nuclei and within QMC itself. 

To summarize, we have calculated the single-particle energies for rj- and tu-mesic nuclei 
using QMC-I. The potentials for the mesons in the nucleus have been calculated self-consistently 
in local density approximation, embedding the MIT bag model rj and u mesons in the nucleus 
described by solving mean-field equations of motion. Although the specific form for the widths 
of the mesons in medium could not be calculated in this model yet, our results suggest that one 
should find 77- and tu-nucleus bound states for a relatively wide range of the in-medium meson 
widths. In the near future, we plan to calculate the in-medium uj width self-consistently in the 
QMC model. 
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Figure 1: Effective masses of the nucleoli, physical rj and u mesons, and those calculated based 
on SU(3) quark model basis (the dotted lines), namely cu = -^(uu + dd) (ideal mixing) and 
r] 8 = -^(uu + dd — 2ss). The two cases for the u meson are almost degenerate. (Normal nuclear 
matter density, p Q , is 0.15 fm~ 3 .) 
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Figure 2: Potentials for the rj and u mesons, (m*(r) — m v ) and (m* s (r) — m WB ), calculated in 
QMC-I for 26 Mg and 208 Pb. For 208 Pb the potentials are also shown for QMC-II. 
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